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Abstract 



We study the doping of a Mott insulator in the presence of quenched frus- 
trating disorder in the magnitude and sign of the magnetic exchange. Two 
quite different doping regimes 6 < 5* and 6 > 6* are found, with 6* ~ J/t 
(J is the characteristic magnitude of the exchange, and t the hopping amph- 
tude). In the high-doping regime, a (Brinkman-Rice) Fermi hquid description 
apphes with a coherence scale of order St. In the low doping regime, local 
magnetic correlations strongly affect the formation of quasiparticles, result- 
ing in a very low coherence scale ep ~ J{5/5*)'^. Fermi liquid behaviour 
does apply below e^, but a "quantum critical regime" e*p < T < J holds, 
in which marginal Fermi liquid behaviour of several physical properties is 
found: NMR relaxation time l/Ti ~ const., resistivity Pdc{T) oc T, optical 
lifetime oc oj / \n.{uj / e*p) together with uj/T scaling of response functions, 
e.g. JJ2qX"iQT^) oc tanh(ci;/2r). In contrast, single- electron properties dis- 
play stronger deviations from Fermi liquid theory in this regime with a ^/uJ 
dependence of the inverse single-particle lifetime and a decay of the 

photoemission intensity. 

On the basis of this model and of various experimental evidence, it is 
argued that the proximity of a quantum critical point separating a glassy 
Mott- Anderson insulator from a metallic ground-state is an important ingre- 
dient in the physics of the normal state of cuprate superconductors. In this 
picture the corresponding quantum critical regime is a "slushy" state of spins 
and holes with slow spin and charge dynamics responsible for the anomalous 
properties of the normal state. This picture may be particularly relevant to 
Zn— doped materials. 

PACS: 71.10.Hf, 71.30.+h, 74.72.-h. Preprint LPTENS 98/24 
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I. INTRODUCTION 



How (and whether) coherent quasiparticles form in a hghtly doped Mott insulator is a 
key question in the physics of strongly correlated electron systems. A satisfactory theoretical 
understanding of this issue has been achieved in the limit where magnetic correlations do 
not play a prominent role, starting with the work of Brinkman and Rice In cuprate 

superconductors however, the undoped phase is an antiferromagnetic insulator with a rather 
large exchange coupling Jaf (on the scale of 100 me V^), so that we have to face the problem 
of the interplay between local coherence and magnetic correlations. 

Furthermore, there is ample experimental evidence that carrier localisation and mag- 
netic frustration also play a crucial role in the low to intermediate doping regime. This is 
particularly clear in the La2-a;Sr^Cu04 compound at concentrations just above x = 0.02 
(the threshold for the disappearance of the antiferromagnetic long-range order), for which 
true spin-glass ordering of the copper moments has been demonstrated at very low temper- 
ature (with Tg ^ 7K for x = 0.04 ||^). Up to which doping concentration does this glassy 
regime persist when superconductivity is suppressed is not known at this point, but carrier 
localization is indeed observed at low temperature up to optimal doping in both the ab and 
c directions when a strong magnetic field is applied It was actually predicted early 

on [0 that hole doping induces strong frustration in the system when the holes become 
localized, replacing locally an antiferromagnetic Cu-Cu bond with an effectively ferromag- 
netic one, with a strength larger than the original Jaf- We observe furthermore that the 
disappearance of antiferromagnetic long-range order is accompanied by the appearance of 
new low-energy spin excitations, of a quite different nature than spin waves, as evidenced 
by inelastic neutron scattering experiments p- [1^. It is important to notice that the com- 
pounds with a glassy ground-state display, at sufficiently high temperature (above the onset 
of localization), the same distinctive transport properties as in samples with higher doping, 



e.g. linear resistivity [§|T0|- It is thus tempting to view these low-energy excitations as the 



source of anomalous scattering in the normal state. 

Anticipating some of the speculations made at the end of this paper, we shall argue that 
these low energy excitations are associated with a novel kind of spin state: the "slushy" 
state associated with the disordering of an insulating (possibly glassy) ground state by 
hole motion, quantum fluctuations and thermal effects. In this picture, many distinctive 
"anomalous" properties of the normal state of the cuprate superconductors are associated 
with the quantum critical regime corresponding to the T = transition at which the insu- 
lating (glassy) ground-state melts into a metaUic (Fermi-liquid) ground state when doping 
is increased. 

In this paper, we shall study a highly simplified model of such a slushy state of spins 
and holes. Our starting point is the work of Sachdev and Ye [jl^, who showed that in the 
large-M limit of the fully-connected random Heisenberg model of SU{M) spins, quantum 
fluctuations are strong enough to overcome the tendency to spin-glass ordering. Instead, a 
gapless spin-liquid state is found down to zero temperature with a large density of low-energy 



spin excitations Remarkably, these excitations are characterized by a local dynamic 

spin susceptibility which has precisely the form advocated by the "marginal Fermi liquid" 
phenomenological description |T^ of the low energy spin excitations in cuprates, namely: 
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XL(^, T = 0) = ^sgn(.;) , ^ = 0) = Ij j 

This model is one of the few cases in which a response function having the marginal Fermi 
liquid form could be derived explicitly (see also The generalization of Eq.(|^) to finite 

temperature will be given in Sec. |III Q (Eq. (0)) and displays u/T scaling. The physical 
mechanism for the gaplessness and the high-density of spin excitations in this model is 
discussed in more detail at the beginning of Sec.|T|. It has to do with the large number 
of transverse components of the spins in the large-M limit. In this respect, it might be a 
reasonable picture for the disordering of the two-dimensional quantum Heisenberg spin-glass 
due to quantum fluctuations and low dimensionality ^1 . 



The main purpose of this paper is to determine whether this marginal Fermi liquid 
spectrum survives the introduction of charge carriers and the associated insulator to metal 
transition. The physics of this problem is dominated by the interplay between two competing 
effects: 

• The formation of coherent metallic quasiparticles, which can be viewed as a binding 
of spin and charge degrees of freedom. In the simplest description of a doped Mott 
insulator with U = oo, coherent quasiparticles form below a scale of order T^q ~ e^Q ~ 
St (where S is the doping and t the hopping amplitude). This is a "naive" estimate of 
the effective Fermi-energy scale, since it ignores any effect coming from the magnetic 
exchange (which will tend to suppress it). 

• The binding of spin degrees of freedom on neighbouring sites into singlet or triplet 
states, and the corresponding slow dynamics of the on-site local moment. This is the 
phenomenon leading to the formation of the spin-liquid state in the undoped phase, 
which involves a scale of order J (the characteristic strength of the exchange). 

It is clear from comparing the scales above that when J is larger than the "naive" 
coherence scale e*pQ, the magnetic exchange prevents the formation of coherent quasiparticles 
at that scale: in other words, e*pQ cannot possibly be the actual quasiparticle coherence scale 
above which free local moments are recovered, since the exchange is still effective at energy 
scales between e^g and J. It is thus expected that the actual coherence scale of the system, 
e*p will be much smaller than e*pQ, and that a new metallic regime in which spin degrees of 
freedom form a spin-liquid like state while charge degrees of freedom are incoherent will be 
found in the intermediate energy and temperature range e*p < uj,T < J. From the above 
estimates, this will be the case at small doping: 6 < 6* ^ J/t, while a direct crossover from 
a coherent metal to an incoherent high temperature state is expected for 6 > 6*. These 
expectations are entirely borne out from our solution of the doped Sachdev-Ye model, as 
evidenced by Fig.|l], which summarizes the main crossovers found in our analysis. 

It should be emphasized that this competition between metallic coherence and magnetic 
exchange is also essential to the physics of heavy fermion compounds |2^. In this context. 



the "naive" coherence scale e*pQ stands for the single-impurity Kondo scale (or rather, any 
estimate of the lattice Kondo scale that ignores RKKY interactions), while J stands for the 
typical strength of the RKKY interaction. For this reason, the results of the present paper 
may also have some relevance, with appropriate changes, to the physics of the disordered 
rare-earth compounds near the quantum critical transition into a spin glass ground state 
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II. THE MODEL 



A. Disordered SU{M) t - J model 

The effect of charge carriers on the Sachdev-Ye spin-hquid phase will be investigated 



by generalizing the model of Ref. to a t — J model, with randomness on the exchange 
couplings Jij between nearest neighbors sites : 

H = ~ ^ tij Pcl^Cja -P + X/ 'Jij'^i ' "^j (2) 
<ij>a <ij> 

In this expression, the SU{2) spin symmetry of the electrons has been enlarged to SU{M) 



25| . Si is the conduction electron spin density on site i and the spin index a runs over 



a = 1, ■ ■ ■ , M. The projection operator P enforces the local constraint: 

. M . . 

2^ cLCia < — (3) 
a ^ 

In this manner the M = 2 case exactly coincides with the standard t — J model with the 
constraint of no double occupancy. 

The exchange couplings are quenched random variables with random sign and magnitude, 
distributed according to a Gaussian distribution with : 

(throughout this paper the bar will denote an average over the disorder). In the following, 
we shall consider this model on a lattice of connectivity z, with a nearest-neighbor hopping 
amplitude normalized as: 



2t 

and we shall analyze the model in the following double limit: 



i) z ^ oo. In this limit of infinite connectivity, a dynamical mean field theory applies 
which reduces the model to the study of a single-site self-consistent problem ||^, as 
detailed in Sec. |lllj|. However this single site model is still a complicated interacting 



problem. For the sake of simplicity, the lattice will be taken to be a Bethe lattice (no 
essential physics is lost in this assumption). 

M — >■ cx), in which the single site problem becomes tractable. In the absence of a 
random exchange, this limit yields the familiar Brinkman-Rice description of a doped 
Mott insulator |IIE5|. 



The scaling in z and M in (^ and (^ are chosen such that this double limit gives non trivial 
results. Alternatively, one could consider (as in [|l^) this model on a fully connected lattice 



of sites, with random hopping amplitudes: t^j = with = 0, ^'fj = 1. This leads 

to precisely the same equations for single particle Green's functions as the z = oo Bethe 
lattice M. 
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We shall use a decomposition of the physical electron operator into a spin-carrying 
fermion / and a slave boson b : cf^ = f^bi. The local constraint (|) becomes : 

E /I/- + = ^ (6) 

With this decomposition the Hamiltonian (0) can be rewritten as: 

2t J 
^ = ~jr7^ ^ (/l^^^j/ja + h.c) + -= E E S^apSjp^ (7) 

and the — 1 components of the SU{M) spin operators Si = {Si)°'^ read : 

1 

SialB = ftafiP ~ TJ'^"/^ E fiafia (8) 



B. Reduction to a single-site problem 

In this section, we explain how the large connectivity limit z ^ oo reduces the problem 
to the study of a single-site model supplemented by a self-consistency condition. First we 
use a path integral representation of the partition function Z and introduce a Lagrange 
multiplier field Ai(r) on each site in order to handle the constraint (]^). We then introduce 
n replicas of the fields (/f , 6^, A" , a = 1, ■ ■ ■ , n) in order to express Z"' and average over the 
disorder. The action associated with reads: 



i,a V^^^^ <ij>,a-^0 

^'^'drdr' H Y. 5L,(^)5;>.(r)5,^,(r')5),,(r') (9) 

<ij> l<a,/3,7,5<A/ l<a,b<n 



2zM Jo Jo 



where the action 5*0 is defined by 

r0 



So[f\f,b\b,X]^ I dr i^b\T)drb{T) + YfMidr - f^)Ur)j 

ri3 ( /If \ 

M I drXir) lYflir)Ur) + b\T)b{r) - - (10) 



Following the "cavity method" (reviewed in 0), a site of the lattice is singled out, and a 
trace is performed over all degrees of freedom at the other sites (concentrating on phases 
without translational symmetry breaking, so that all sites are equivalent). In the z oo 
limit, this can be performed explicitly, and the problem reduces to a single-site effective 
action which reads: 

Seff = E ^o[/^^ &^ ^1 - ^ E /7' drdr' ^:,(rX,^(r - T')S'^sir') 

+ Q)'E drdr' fl\T)b\T)CZ{r - rO&^"(rO/a"(rO (11) 
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This effective action is supplemented by a self-consistency condition which constrains C{t — 
t') and R{t — t') to coincide with the local electron Green's function and spin correlation 
function respectively: 

CZiry) = -{Tct{r)ct{r'))^ = -{T{rJ^){T){fi%'^){T')) ^^^^ 

K^sir, r') = {Sl,iT)S':^sir'))^ = (5^,(r)5!;,(r'))^^^^ (12) 

In each of these equations, the last equality expresses the fact that local correlation functions 
can be calculated using the single-site action S'e// itself. The limit n ^ must eventually 
be taken in these equations. 



C. Saddle-point equations in the large-M limit and slave-boson condensation 

We shall study the above self-consistent single-site problem in the large M limit, focusing 
on the paramagnetic phase of the model. In this case, all the above correlators become replica 
diagonal (C'^'^ = C, D""^ = D5ab)- 

Furthermore, we shall look for solutions in which the slave boson undergoes a Bose 
condensation. (Solutions with an uncondensed boson when the bosons carry an additional 
channel index have been investigated by Horbach and Ruckenstein p4l). The solutions 



considered here can be found as a saddle-point of Seff by setting: 6(r) = yM/20(r) and 
looking for solutions in which both 0(r) and the Lagrange multiplier A(r) become static at 
the saddle-point: 

hspir) = 5 ^Asp(r) = Aq (13) 

From the constraint Eq.(^, the total number of electrons will be related to 6 through: 
J2a < flfa >= ~ ^'^ ^^^^ ^ measures the number of holes doped into the system. 

The saddle point equations then reduce to a non-linear integral equation for the fermion 
Green's function — < T/a(T)/j(r') >= dapG f{T — t'), which reads (with : Un = {2n+l)7T/ (3 
the Matsubara frequencies): 

Cj^iiujn) = iujn + - Xo - itSfGfiiujn) - S/(iu;„) (14a) 



-j'G){r)Gf{-r) (14b) 



and to the following relations, which determine the Lagrange multiplier Aq and the chemical 
potential /i for given values of the doping 5 and the temperature (given the couplings J and 
t): 

Gf{T = 0") = ^—^ (14c) 



XoV6 = -2f6-^ [^dTGf{T)Gf{-T) (14d) 

<^ 
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The derivation of these saddle point equations from Sejf is detailed in Appendix 0. 

The local spin-spin correlation function is directly related to in the M — oo limit, 

as: 

Ri^) ^^T.< S^o.p{T)S,p^{0) >= -Gf{r)Gf{-r) (15) 

a/3 

In the following, we shall often consider the spectral functions associated with the single- 
particle Green's function and the local spin-spin correlation: 

Pf{u) = --lmGf{u; + zO^) , xL(^) = Im + zO+) (16) 



III. PHYSICAL PROPERTIES OF THE METALLIC STATE 

In this section, we study the nature of the metallic state as a function of the doping level 

5. 

Let us first recall some of the properties of the spin-liquid insulating state found for 
5 = 0, as obtained by Sachdev and Ye |]T^. In this case, our equations ( |14a| - |14cD coincide 



with those of Ref. |T^. Note that Eg .( 114011) decouples, being automatically satisfied at 5 = 0, 
and that particle-hole symmetry imposes /i — Aq = Sj(zO"'") = 0. A low-frequency analysis 
of the integral equation reveals that the T = Green's function and spectral density have a 
I/a/cJ singularity for |ci;| — >■ 0. More precisely ||2^, in the complex frequency plane as 2; — 0: 



G,(.)=(^)'''^ + --- , lm.>0 (17) 



This yields the following behaviour of the local dynamical susceptibility for u ^ 0: 



4J 



sgn(tu) + ■ ■ • (18) 



Fig.|| displays a numerical calculation of Pj{oj) and Xiocl^) ^^ro doping (in agreement 
with the one in Ref. [T^). These results display the above low-frequency behaviour (but we 
note that significant corrections to (|18D are already sizeable at rather low values of uj / J .) 

Hence the insulator at 5 = is a gapless quantum paramagnet (spin liquid), with a 
rather large density of low-energy spin excitations. Remarkably, (|T8|) is of the same form 
than the "marginal Fermi liquid" susceptibility proposed on phenomenological grounds by 



Varma et al. |19| for the normal state of the cuprate superconductors. In the present context, 
the physical nature of these low-energy excitations is intimately connected to the fact that 
the exchange couplings Jij are random in sign. In constructing the ground-state of the 
insulator, let us imagine that we first try to satisfy the bonds with the larger exchange 
constants. When such a bond is antiferromagnetic, the two spins connected by it will form 
a non-degenerate singlet. For a ferromagnetic bond however, the two spins will pair into a 
state of maximal possible spin (the generalization to SU{M) of a triplet state). This state 
has a degeneracy, which actually becomes very large (exponential in M) as M becomes large. 
Continuing the process in order to accommodate bonds with smaller strengths will tend to 
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remove part of this degeneracy , but leaves behind a very large density of low-energy spin 
excitations. These effects are clearly favored by the (fermionic) large-M limit considered 
here, because of the high degeneracies of the "triplet" state and because the strength of 
quantum fluctuations in this limit precludes the appearance of long-range order (e.g. spin- 
glass) which would remove degeneracies in a different manner. We believe however that this 



physics is not an artefact of the large— M limit. Indeed, preliminary theoretical studies p8 
suggest that the local spin correlations near the quantum critical point associated with the 
T = transition into a metallic spin-glass phase could be similar to Eq. (p!8|) , with a related 
physics. 

Finally, we note that the single-site action to which the model reduces at zero doping 
{i.e Eq. (|llD with t = 0) has some similarities with the multichannel Kondo effect in the 
overscreened case. In the present context however, the "bath" seen by the spin is not due to 
an electronic conduction band, but generated by all the other spins in the lattice. The spin 
correlations of both the bath and the spin adjust to the self-consistent long-time behaviour: 
< S{0)S{r) >~ R{t) ~ 1/r similar to that of the SU{M) Kondo model with K = M 
channels BO] (2-channel model in the SU{2) case). 



A. Low-frequency analysis and the Fermi liquid coherence scale 

The first question we would like to address is whether the "marginal Fermi liquid" spin 
dynamics survives the introduction of charge carriers. As we shall demonstrate, this depends 
on the temperature range considered (Fig.|l|). At low temperature, below some -possibly very 
low- coherence scale e*p, it turns out that a Fermi liquid is recovered. 

This is easily seen from a low-frequency analysis of the integral equation for at 
zero-temperature. At zero doping, the Green's function and self-energy behave at low- 



frequency as: Gf{uj) ~ 1/yJuj , Sj(c<j) ~ yJuj. When inserted in Eq. ( |14a| ) , this controls 
the leading low-frequency behaviour of both the r.h.s and l.h.s of the equation taken at 
(5 = 0, which match each other. However, for 5 7^ 0, the term {t5YG f{uj) would introduce a 
'^1 \/lo singularity and prevent this matching from taking place: this indicates that the low- 
frequency behaviour of the zero-temperature Green's function for arbitrary small doping is 
no longer I/a/cJ- In this respect, an infinitesimal doping is a singular perturbation of the 
above equations. This observation directly yields an estimate of the coherence scale e*p such 
that Gf{uj) ~ is recovered for e*p < uj J. Indeed, the term {t6)'^G f{uj) becomes 

comparable to ^f{uj) in this regime (thus providing a cut-off to the singular behaviour) when 
u ~ {5tY / J. Hence, in the low-doping regime: 

4 = ^ , (5 « n (19) 

where 5* will be precised below. (In the following we shall take Eq. (p^ ) as defining e*p in the 
low-doping regime, with no additional prefactors). 

In the high-doping regime on the other hand (or when t 3> J), one should consider 
first the limit of a vanishing magnetic exchange J = 0. In this limit, the usual slave-boson 
(large M) description of a doped Mott insulator is recovered Setting J = E/ = in the 
equations above yields a semi-circular spectral density: 
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where D is given by : 



pf" - \o (20) 



(21) 




The original bandwith 4t of the non-interacting case has been reduced by a factor 5, and 
the usual Brinkman-Rice result for the coherence scale is recovered : 

e*p = t5 , (5 > 5*) (22) 

Turning on J as a perturbation from this starting point does not affect the leading 
low-frequency behaviour of the self-energy, but does lead to a scattering rate ImS/ oc 
uj'^J'^/{5tY + ■ ■ ■ characteristic of a Fermi liquid (in contrast the J = model has infinite 
quasiparticle lifetime in the large-M limit). From Eqs.( p^JT9| ), it is clear that when the mag- 
netic scattering is strong (J ^ t), regime ([19D always applies, while for weaker scattering 
{J < t) a crossover between the two regimes is found at a characteristic doping: 

6* ~ min(-, 1) (23) 

t 

We thus observe that below some characteristic doping the low-energy coherence scale is 
strongly affected by the magnetic scattering. When the exchange is large or for doping 
smaller that 6* ~ J/t, the actual coherence scale e*p is much smaller than the "naive" 
coherence scale e*pQ (which holds in the absence of magnetic correlations). Here we find 
e^Q ~ 6t and e*p/e*pQ ~ 6/6*. This is one of the crucial physical conclusions of this paper. 

A numerical solution of the saddle-point equations provide a clear evidence for these two 
regimes. The numerical procedure that we have used is explained in Appendix |^. Fig. (^ 
displays the T = spectral function Pf{uj) for three values of J. When J is very small, the 
spectral function is very close to the semi-circular shape (pO]), while for a larger J the ^ 
divergence is observed over a large frequency range e*p < uj < J but is cutoff for c<j < so 
that p/(0) is finite. Anticipating on the results of Sec. |111 B| , we observe that the value of 



pf{uj = 0, T = 0) is actually independent of J as a consequence of the Luttinger theorem. 
Indeed, the following relation can be established at zero-temperature: 

/i- Ao-Sj(iO+) ^5/io(5) as T^O (24) 

where po{6) is the non-interacting value of the chemical potential for the tight binding model 
on the z = oo Bethe lattice. This implies: p/(0,T = 0) = l/(7rt5) for all values of J. 

At very small doping, a scaling analysis of the saddle-point equations can be performed 
in order to characterize more precisely the crossover between the low-frequency and high 
frequency regimes at T = 0. As we now show, the spectral function (and the Green's 
function itself) obeys a scaling form: 

P/M = ^0/(j) for a;«J,t ^ « ^* = J (25) 
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In order to derive the integral equation satisfied by the scahng function 0/, we rewrite 
Eq.(|l4i) at T = (using (H) as : 



G]\co) =uj + 5fXoiS) - itS)^Gfiuj) - (Sj(cu) - S^(0)) (26) 

At low doping and low- frequency G/ is of order 1/6 and /iq is of order 6. Hence, rescaling 
frequencies by the coherence scale e*p = {Sty/ J, we see that the first two terms in the r.h.s of 
(p6|) can be neglected. Analytically continuing to realtime t and frequency u, and denoting 
by the real- frequency T = Green's function (with the usual Feynman prescription), 
we define a scaling function gj associated with Gj by the Hilbert transform: 



oo 



9f{u^)=r'de '^f.^ (27) 
■' J-oo uj — e + I sgntu 

We finally obtain from ( p6D an integral equation satisfied hj Qf (and thus by (pf) which no 
longer contains dimensional parameters: 

{gf{uj))-' = -gf{uj)-{a^{uj)-a^{0)) 

^^(t) = {gf{t)rgf{-t) (28) 

(as explained in Appendix. |C|, a sign change occurs in the expression of the self-energy at 
T = 0) 

The universal scaling function (pf can be obtained by solving numerically Eq. (28), and 
the result is displayed in Fig. ^. The asymptotic behaviours of (f)f{uj) for large and small 
UJ = uj/e*p can be obtained analytically and read: 

(pfip) = ^ — Ciuj'^ + ■ ■ ■ for ttJ — > 

(Pf{uj) = ^ + . . . for u; -> +CX) (29) 



where ci and C2 are two constants. The low-frequency behaviour reflects the Fermi-liquid 
nature of the low-energy excitation spectrum, while the \/\JZj behaviour characteristic of 
the undoped spin-liquid is recovered for uj > e*p. 



B. Single electron properties at T = : quasiparticle residue, effective mass, 
Luttinger theorem and photoemission 

In this section, we focus on the one-particle Green's function for the physical electron, 
which is related to that of the auxiliary fermion by: 

G,(k, luj^) = -< Tc^^cl^ >=-< Tblf^jlj^ >= —Gf{k, z^„) (30) 



hence: 

M5 



GciKiuJn) ^ = icun + /i - Ao - S/(2t^„) - (5e]^ (31) 
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In this expression, ej^. stands for the one-particle energies of a non-interacting tight-binding 
model on the Bethe lattice with hopping t/y^ between nearest neighbour sites |3l|]. The 
distribution of these single-particle energies is the semi-circular density of states D{e) defined 
in(P). 

From the large- frequency behaviour of Eg. (|31|) , we see that the physical electron spectral 
density in the M ^ oo limit is normalized as pc = M6/2 (in contrast pj = 1). This 
is expected from the fact that the constraint on the Hilbert space yields a normalisation 
(for arbitrary M) J^^ Pc =< {c,c^} >= M6/2 + (1 - S)/2 (note that this yields (1 + 6)/2 
for M = 2, as expected for the U = oo Hubbard model). 

Since our normalisation of the hopping is Uj = 2t/{My/z), the non-interacting conduction 
electron Green's function reads: 



2 

Gc{k,iuJn)]ree = i^n + P- J^e^. (32) 



Thus, the physical electron self-energy reads: 



2 

Ec(iu;„) =iujn + fJ'- jj^ {i^n + p - Xq- T.f{iuJn)) (33) 

We observe that it depends solely on frequency, as is generally the case in the limit of large 
dimensionality f^. 

We first consider the location of the Fermi surface for both the non-interacting and 
interacting problems, i.e look for the poles of the electron Green's function. In the non- 
interacting case, we relate the chemical potential at T = to the number of particles 
< Ha >= (1 - 5)/2 and find: 

2 

/i/ree = ^/io('5) (34) 

where the function po{S) is defined by the relation: 

m(<5) I- S 

deD{e) = ^ (35) 

-oo Z, 

Hence the non-interacting Fermi surface corresponding to a doping 5 is defined by ej^ = 
Pq{5). In the interacting case, we see from Eq.(^) that the Fermi surface is located at 
~ = 0) ~ ^o(^ = 0) ~ ^/(^ = 0,T = 0)) /5. In the absence of magnetic scattering 
(J = 0), it can be shown by an explicit calculation [0] from the saddle-point equations that 
the r.h.s of this equation is just po{5) and thus that the Fermi surface is unchanged in the 
presence of the constraint. When J 7^ 0, such an explicit calculation is not possible, since 
the saddle point equations are coupled non-linear integral equations. However, a proof of 
Luttinger theorem can still be given using the fact that a Luttinger-Ward functional exists 
for this problem and is known in explicit form in the large-M limit, as detailed in Appendix 
0. The conclusion of this analysis is that the volume of the Fermi surface corresponds to 
(1 — (5)/2 particles per spin fiavour and that the zero temperature, zero frequency self energy 
must obey : 

p{T = 0) - Ao(T = 0) - Sj(cu = 0, T = 0) = 5pq{5) (36) 
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We now consider the weight and dispersion of the quasiparticles, that can be read off 



from Eqs.(31,33) by expanding around the Fermi surface. We define a renormahsation factor 



for the auxihary fermions as: 



1 



(37) 

1^=0 



so that the physical electron quasiparticle residue reads: 

Zc = Y ^^f (38) 

From the low-frequency analysis of the preceding section and the corresponding estimates 
of the coherence scale, we expect Zc to be of order e*p/t and thus : 



2 



Z^r^id^ r^^ (5 < 5*) , Z,r^6 (5 > 6*) (39) 

In FigJ^, we display the result of a numerical calculation of Zc as a function of doping, for 
three values of J/ 1. These results entirely confirm the above expectations. We have checked 
that at small doping Zc/S* scales proportionally to {5/5*Y with a universal prefactor. 

From Eq. (pT]) , we see that the quasiparticles have a dispersion characterized by an effec- 
tive hopping teff/t = SZf {m* /m = l/Z^ oc l/{6Zf)). Hence the effective mass diverges as 
the Mott insulator is reached (as 1/5^). The reason for this divergence is the large (extensive 
27| ) entropy of the insulating spin-liquid ground-state. This entropy must be released at a 



temperature of the order of the coherence scale e*p in the doped system. Hence, integrating 
the specific heat ratio C /T = j between T = and T = e*p leads to 76^ ~ 1, which is the 
result found above. This divergence of 7 as 5 —> is clearly an artefact of the large-M and 
large-c? limit. The residual ground-state entropy of the spin-liquid phase should not survive a 
more realistic treatment of this phase (whether this happens while preserving x"{^) ~ const, 
in this phase is an open problem at this moment). Furthermore, our model does not include 
a uniform antiferro magnetic exchange constant superimposed on the random part. Includ- 
ing this coupling will help locking the spins into singlets and cutoff the divergence of the 
effective mass (for a large-M treatment of this point, see e.g Ref. 0). 

Finally, we discuss the shape of the conduction electron spectral density p(ejj.,co') for 
a fixed value of the energy ej^. as a function of frequency, as relevant for photoemission 
experiments : 



. ^ ^ MS s^M 

Pc(ek' ^) = (ek- ^) = ^2 (40) 



Numerical results for this quantity are displayed on Fig. ^. This function is peaked at a 
frequency ujpeak — -^c(ek~^k ) ^i^h a height of order l/ci^pg^^ (at T = 0). Moving away from 
this quasiparticle peak, pc{€-^^,u) has the characteristic l/cu^ decay of a Fermi liquid only 
in the limited frequency range jcUpeafcl < 1"^! < ^f, followed (for 6 < 6*) hj a much slower 
1/y/uj tail corresponding to the spin liquid regime in the frequency range e*p < \uj\ < J. 
(We note that this non- Fermi liquid tail is absent in the high-doping regime). These two 
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regimes are clearly apparent on Fig. p. If the resolution of a photoemission experiment is 
not significantly smaller than e*p, the peak will be smeared into a broad feature, and the 
measured signal will be dominated by the slowly decaying tail. Furthermore, as shown in the 
next section, temperature has a large effect on the peak, whose height decreases as l/VT 
in the temperature range e*p <T < J. 



C. Finite-temperature crossovers 

The metal-insulator transition at T = as 5 — is a quantum critical point. The 
associated crossover regimes at finite temperature can be easily deduced by comparing the 
coherence scale to the magnetic exchange J and to the temperature. This analysis yields 
three regimes, as depicted on Fig.|l]: 

• For T < e*p, the doped holes form a Fermi liquid. The low-energy degrees of freedom 
are the fermionic quasiparticles described by the auxiliary fermions /q,, which behave in 
a coherent manner since their inverse lifetime vanishes at low-frequency as Im S/ oc a;^ 
in this regime. 

• At low doping 5 < 5*, an intermediate temperature regime exists, defined by: e*p < 
T < J. In this regime, coherent quasiparticles no longer exist (as shown below, 
ImSj oc \/uj), but the spin degrees of freedom are not free local moments since the 
temperature is smaller than the magnetic exchange. Hence, the spins behave in this 
regime as in a spin liquid, with a marginal Fermi liquid form for the local spin response 
function. As shown below, this regime corresponds to the so-called "quantum-critical" 
regime associated with the quantum critical point at T = 5 = 0. In this regime, the 
low-energy scale e*p drops out from response functions, which obey universal scaling 
properties as a function of the ratio uj/T. 

• Finally, a high-temperature regime applies, defined by T > J (for 6 < 6*) ot T > 6t = 
e*p (for 6 > 6*) in which both spin and charge are incoherent and essentially free. We 
note that if J < t and the doping is larger than S* ^ J/ 1, the system goes directly from 
a Fermi liquid to this high-temperature regime as temperature is increased, without 
an intermediate marginal Fermi liquid regime. 

This qualitative analysis can be established on firmer grounds by generalizing the low- 
doping scaling analysis of Sec. |lll A| to finite temperature. Assuming that the coherence scale 



is small as compared to both J and the hopping t (i.e that S < S*), and that uj,T <^ J,t, 



the spectral function takes the following scaling form, generalizing Eq.(|25D: 




P/(^,T) = >W^,^ (41) 



In the following, uJ and T stand for uj/e*p and T/e*p, respectively. Eq. (|4l|) yields for the 
Green's function Gf{uj) = j^QfioJ^T) (with $/ = — (^j/vr). We assume that the self energy 
also scales as Sj(co') = (t5)a"^{jjj,T). From the saddle-point equation (|^) we deduce the 
following equations : 
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Im^?7'(^,T) = -g'}(p,T) - a"f{uJ,T) (42) 

/oo roo 
/ dXidX2 ^ f{Xi,T)^ f{x2,T)^ f{Xi + X2 -uJ,T) 
-oo J — oo 

In this expression, np and ub are the Fermi and Bose factor : np^Biv) = l/(e^ ± 1). With 
a Kramers- Kronig transformation one can deduce g'j- and o"j(cJ, T) — a'j{0,T) from and 
(Tj. From the equation for RegJ^, we have that /i — Aq — Sj(ct; = 0,T) ~ f{T/e*p) where 
/ is some scahng function which can in principle be calculated from (|42|). The function / 
vanishes at small argument (/(O) = 0) due to Luttinger's theorem, and at large argument 
(/(+oo) = 0) due to the auxiliary fermion particle hole symmetry of the undoped model. 

We now discuss the solution of this scaled integral equation, and the form taken by 
in the various regimes described above. 

i) Fermi liquid regime : T <^ e^. At zero temperature, the scaling function $j reduces 
to that in Eq.dl): 

$/(cJ,T = 0) = (44) 

We can also consider the limit of low-frequency and temperature uj,T <^ e*p but with an 
arbitrary ratio u/T. In this limit, the self-energy term is negligible altogether in Eq.(^), 
and one gets simply gj = —1 i.e: 

$/(cJ< 1,T< 1) ^ - (45) 

TT 

Note that the r.h.s could a priori be a function of the ratio u/T, but is actually a constant 
(as is generically the case in a Fermi liquid). From this, we can deduce a scaling form of the 
scattering rate in the same regime. Indeed ( ^5]) corresponds to the imaginary time Green's 
function: 

Gfir)^-^-^^ 1/6^ «r,/3-r (46) 

Trtdsmnr/fj 

Hence, in this limit, the self-energy takes the form: 

^^("^ ^ - (^ (^i^) (^^^ 
which can be Fourier transformed to yield: 

ImS;(^ « e^,T « e^) = -^^i^' + ^'T') (48) 

ii) Spin-liquid regime : T ^ e*p. In this quantum critical regime the energy scale drops 
out from the problem and the spectral density and response functions become functions of 
the ratio oo/T only. Indeed, the scaling function $/(a;, T) takes the form (pf{uJ/T)/\/T in 
the limit T ^ 1. In order to find ipj in explicit form, we divide both side of ( ^21) by Vt and 
take the limit T — > oo, uJ/T fixed. Then the first term of the r.h.s vanishes and we are left 
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with a scaled equation for (pf in which all dependence on e*p has disappeared. Remarkably, 
this integral equation can be solved in closed form and yields : 



Pf{uj,T) 



iff 



27rl- 



Or 



: cosh 







■ ^ \ 


(2t) 




h I 






2nTj 



(49) 



Some details are provided in Appendix 0. This scaling function describes how the 
singularity associated with the low-energy excitations of the spin-liquid is cutoff (by the 
temperature) at frequencies u; < T so that the spectral density is of order l/y/ JT at u = 0. 
(Note that if the limit T — > cxd is performed while keeping uj fixed in (0), the same result 
is obtained as when the limit is taken with uJ = 0. Hence there is no additional crossover in 
the frequency dependence of the response functions below u = e*p in this regime). Eq.f^ 
corresponds to the following scaling form for the imaginary time Green's function: 



Gfir) ~ 



1 



1/2 



2J'Ki \simTT/f3^ 



(50) 



Remarkably, ([50| ) has the form which would hold in a model having conformal invariance, 
for example a quantum impurity model of a spin interacting with a structureless bath of 
conduction electrons. In that case, a conformal mapping from the T = half plane r > to 
the finite-temperature strip < r < P can be used to show |^ that if the Green's function 
decays as at T = 0, then it takes a scaling form given by (|50| ) at finite temperature 

(lower than a high-energy cutoff). In the present case, the original model is an infinite 
connectivity lattice model which does not a priori satisfy conformal invariance. It does map 
onto a single-site quantum impurity model, but with an additional self-consistency condition. 
This means that the effective bath for the local spin is given by the local spin-spin correlator 
itself, and thus does have non-trivial structure at low-energy. However, this structure appears 
only as a subdominant correction to the leading low-frequency behaviour x'l'od^) ~ const. 
For this reason, our effective single-site model does obey conformal invariance properties 
in the low-energy limit, which explains the result above. This remark actually applies in a 
broader context than the specific model considered here, as will be discussed in more details 
elsewhere. 

Let us also consider the scattering rate in this regime, which is obtained by Fourier 
transforming the imaginary time self-energy: Sj(r) ~ — a/J ( sinir/fi ) ^ /(47r)i which yields: 



S}(cj) ~ -vr-tv^cosh 



V2T 



3 . UJ \ 

4 + ^2^j 



(51) 



We have calculated numerically the real-frequency, finite-temperature Green's function 
by following the method described in Appendix |C[ On Fig. ^ we display results for the 
spectral density for various temperatures for J/t = 0.3 at a doping of 5 = 0.04 < 5*. 
These values correspond to a low energy coherence scale e*p/J = (6/6*)'^ ^ 1.810"^. The 
crossover from the Fermi-liquid regime at low temperature into the quantum critical regime 
at intermediate temperatures is clearly visible (in particular, the peak height can be checked 
to decrease as l/VT). Note also that Pf{Lo) remains approximately centered at ^ until 
T c:^ J and shifts rapidly away from uj = for T > J. In the inset, we also display the 
thermal scaling function associated with pf, Eq. (^9]) . 
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D. Local spin dynamics 

In this section, we describe the behaviour of the local spin dynamics in the various 
temperature regimes. In the large-M limit, the local spin correlation function is given by: 

xUr) = -Gf{r)Gf{-r) 

/+00 
dupf{u)pf{u - uj) {npiiy -uj)- npiu)) (52) 
-oo 

In Fig. ^, we display x'Lci^) for various temperatures and the same choice of parameters 
as in Fig. |^. In the low doping regime, e*p <^ J,t, 6 < 6* , Xioc obeys a scaling form which 
follows from the convolution of fffll) : 



XL(-.T)4* (53) 

Let us discuss the limiting forms of this expression as T ^ or T ^ e^. 

i) At zero temperature, Xiod^) ^ shape which resembles the undoped spin- liquid 
case (Fig. |^) for frequencies uj > e*p. At lower frequency, the Fermi liquid behaviour 
Xioci^) oc a; is recovered. This results in a peak with a height of order 1/J at T = 0. 
This crossover can be described by a scaling function : 

XL(^« J,T = O) = ^0,(^^^ (54) 

where ~ ^xi^^V ~ 0) '^^^ be obtained by convoluting 0/ with itself resulting in 

the asymptotic behaviours : 

xL(u;,T = 0)^^ ^«6^ 

2^ < < J (55) 



3 



This can be used to estimate the behaviour of the static local susceptibility at low 
doping : x'lod^ = 0) = / duxiod^)/^- ^bis integral, the region e*p < u < J 
(corresponding to spin liquid excitations) gives the dominant contribution, leading to 
the logarithmic behaviour for 5^5*: 

xLi^ = 0) ^ i In 1^ (56) 

In contrast, as detailed in Appendix ^ the uniform static susceptibility x = x'(q = 
0,1^7 = 0) is a constant of order 1/J, with no divergence at small doping. 

ii) In the spin-liquid regime T ^ e*p, x'loc becomes a function of uj/T. The corre- 
sponding scaling function is remarkably simple : from Eq. (|50D we have x«oc(t) oc 
7r/(/5sin(7rr//3)) which yields : 



x'LMT) = {^tanh^ (57) 

This behaves exactly as the spin response function postulated in the marginal Fermi 
liquid phenomenology |jl9| {uj/T for uj << T, const, for uj > T). 
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We finally use these results to compute the temperature dependence of the NMR relax- 
ation rate : 

TiT UJ 

Expanding the scaling form (^) to linear order in u (and noting that ^^{0,y) = because 
x" is odd), we get for T <^ J: 

(with ipi^y) = ydx^^{x = 0,y)). In Fig. ||, we plot this universal scaling function. We have 
also checked the data collapse of our numerical results on this function. Limiting forms are 
easily obtained from Eq. (|55D and (0) : 

i) T <^ e^: ■ip{y ^ 1) ~ y/vr Hence ^ ~ ^^^^ . We find a Korringa law (as expected 
from a Fermi liquid) but with a very strong doping dependence. We also note that in 
contrast to a non interacting Fermi gas, l/(TiT) oc l/(t6y, Xioc oc l/J\n{6*/6) and 
X oc 1/ J obey quite different behaviour as a function of doping. In particular the 
so-called "Korringa ratio" l/{TiTx^) — {S*/S)^ is very large at low doping. 

ii) J > T ^ e*p: ip{y ^ 1) ^ hence ^ ~ 47 = const, as expected in a marginal 
Fermi liquid. We note that 1/Ti is doping-independent in this quantum-critical regime. 
This is because the scale e*p no longer appears explicitly. 



(58) 

aj=n 



E. Transport and frequency-dependent conductivity 

In the limit of large connectivity, the current-current correlation function has no vertex 
corrections, due to the odd parity of the current (See e.g. [Q). Hence the frequency dependent 
conductivity is given by : 

Rea{u) = e r deDie) T dup,{e, z/)p,(e, u + ^ fF{i^)-np{i^ + uj) ^^^^ 

J —00 J —00 (jj 

where pc{e, lj) is the single-electron spectral density defined in Eq. (|40|) . This expression yields 
the conductivity in units of e^/{ha'^~'^) where a is the lattice spacing and some numerical 
prefactors have been dropped (we shall also ignore the prefactor M in p^)- 



1. Resistivity 

We first discuss the behaviour of the dc- conductivity : 



/oo /"OO rjij' 

deD{e) / plie,Tx) (61) 

^00 -^-oo 4cosh (^f 1 

i) In the Fermi liquid regime T e*p, we have from the behaviour (|^) of the scattering 
rate : -ImS/(cj,T) oc J'^{uj^ + Ti'^T^)/{5tf and + p - ReS/(u;,T) = uj/Zf + const. T. 
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Making the change of variables e = Tu, we see that the integral over u in adc/T diverges as 
Hence, we find in this regime the expected Fermi-liquid behaviour of the resistivity 
Pdc = l/o"dc: 




p,,(r) (X - T « (62) 



a) In the spin liquid regime e*p <^T <^ J at low doping, —T,'j(uj) is of order V JT (times 
a scaling function of uj/T). This must be compared to 6e ~ 5t in the denominator of pd^, 00). 
Since T 3> e*p, we see that ReSj always dominates over 6e which can thus be neglected. 
Hence one can replace pc{e,uj) by the local spectral function dpfiu). In other words the 
limit (5 — i> must be taken before the low temperature limit in this quantum critical regime. 
Using the thermal scaling function Eq. (^Q]), we obtain : 



A2 ^00 flj. 



The integral can calculated exphcitly using dx F f ^ + ix) = tt^ ( |^, Eq. 6.412), we 
finally find : 

p,e(T) = 16v^ ^ e*p<^T<^J (64) 

Hence the resistivity turns out to have a linear behaviour as a function of temperature in 
the spin-liquid regime, again as in the marginal Fermi liquid phenomenology. This is rather 
remarkable in view of the fact that the single-particle scattering rate behaves as \/T in 
this regime. As further discussed in the conclusion, this is characteristic of a regime of 
incoherent transport in which the transport scattering rate cannot be naively related to the 
single-particle lifetime. Furthermore, we note that the y/u behaviour of the self-energy is a 
crucial ingredient in producing a T-linear resistivity. With a different power law (cu"), the 
resistivity would behave as T^" in this incoherent regime. 

The crossover from to T in the resistivity can be captured in a more precise manner 
in a universal scaling function : 



p.c(T) = ^J-] (65) 




We have determined numerically the function ipp, which is depicted in Fig. |T0|. We observe 
that it is linear over a wide temperature range (with a slope in agreement with 



2. Optical conductivity 

We now turn to the analysis of the frequency- dependent conductivity, 
i) In the Fermi liquid regime, the conductivity takes the form at T = : 

(t(u;) = D6{lj) + aregi^) (66) 
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where D is the weight of the Drude peak and areg{^j~>) — > const, as c<j ^ 0. The Drude peak 
is easier to capture by a finite temperature analysis : the delta function is regularised by 
T in the form T'^ jiixP' + T^). Performing a low-temperature, low-frequency analysis of ( |60[ ) 
leads to the estimation D oc t'^D{fiQ)Zf6 oc 5^ at small doping. A closed formula can be 
given for Kearegi^) as (truncated) convolution of the scaling function (pf. A low-frequency 
analysis then shows that Rea{uj <^ e^) = const. , while Rea{uj ^ e^) ~ e*p/uj. 
a) In the regime e*p < T < J, we have from (pDf) and the scaling form : 



Rea(^) = —ip„ - 

UJ \I 

/+00 ijqr 
, , „ ^/ {^V) ((1 + x)y) [/ {xy) - / ((1 + x)y)] (67) 
^|x(l + x)| 

where /(x) = l/(e^ + 1). From (|67D, (y9o-(+C)o) = const, and thus we have in this spin- liquid 
regime : 

e* 

Re a{uj) oc — T <^ lu <^ J 

UJ 

Rea(cu)oc^ w<T (68) 

Moreover, using the Kramers Kronig relation Imcr(ci;) = / (iti;'Recr(ti;')/(a; — u') we find in 
the same regime for uj > T : 



Ima(cu) oc ^In f ^) (69) 



Hence defining an optical scattering rate from an effective Drude formula : T^p^{uJ) = 
uj'Rea{ijj)/ \ra.a{uj) we find T~p^{|JJ) ~ uj/ ln(u;/e^). 

We have also calculated Recr(a;) numerically, following the method explained in Ap- 
pendix |C[ Numerical results are displayed for various temperatures on Fig.|TT] and are in 
agreement with the previous analysis. 

IV. CONCLUSION AND DISCUSSION 
A. Summary 

In this paper, we have solved a model of a doped spin-fiuid with strong frustration on the 
exchange constants Jij. The undoped model is an SU{M) quantum Heisenberg model with 
random exchange, previously studied by Sachdev and Ye |]T^ in the limit of large-M and 
infinite connectivity. These authors found that, in this limit, quantum fiuctuations are so 
strong that no spin glass phase forms |]TB|. Instead, a gapless spin liquid is found with local 



spin dynamics identical to the marginal Fermi liquid phenomenology |T^. We generalised 



this result to finite temperature and found that the local spin response function displays 
uj/T scaling: Jx" {uj,T)ioc oc tanhu;/2T (for uj,T < J). Doping this Mott insulating phase 
with holes, we found that a characteristic doping 6* ~ J/t appears separating two quite 
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different doping regimes. In the high-doping regime 5 > S*, magnetic effects are weak and a 
Brinkman-Rice Fermi-hquid description is vahd, with a rather large coherence scale of order 
6t. In the low doping regime however, the interplay between local coherence and magnetic 
effects gives rise to a coherence scale e*p = J{5/5*Y^ which can be very low. At low temper- 
ature T < e*p, Fermi liquid behaviour is recovered, but an incoherent regime is found in a 
rather wide regime of temperature e*p <T < J 'm which physical properties strongly deviate 
from Fermi liquid theory. This regime corresponds to the "quantum critical regime" asso- 
ciated with the metal-insulator transition which in this model happens dX 5c = T = {). We 
found that both transport properties and response functions in this incoherent regime behave 
as in the marginal Fermi hquid phenomenology, namely pdc oc T, Topt{oj)^^ oc cu/ ln(co'/e^), 
1/Ti oc const., and J'x'l^^uo^T) oc tanhct;/2T. Remarkably, single-particle properties deviate 
much more strongly from Fermi liquid theory, with a single-particle scattering rate behaving 
as ImE oc ^/uj (or -\/T), in contrast to the ImS oc behaviour postulated in the Marginal 
Fermi Liquid phenomenology. 

These behaviour result from the solution of the large-M saddle point equations, which 
also yields explicit expressions for the scaling functions of uj/e*p and uj/T describing the 
crossover of the various physical quantities between the Fermi liquid and the non-Fermi 
liquid regime. We also note that in the large-M limit, response functions can be calculated 
from the interacting single-particle Green's function. Hence the behaviour of S oc ^/uJ 
and of ImXioc oc const, are intimately related. In contrast, in the marginal Fermi liquid 
phenomenology, the behaviour of Imx is related to a priori unknown higher order vertex 
functions. In this sense, the present model yields a solution to the problem of internal 
consistency of the Marginal Fermi Liquid Ansatz, resulting in a more singular form of the 
single-particle Green's function. 

We also note that numerical studies of the doped two-dimensional t — J model with 
uniform antiferromagnetic J by Imada and coworkers have some intriguing similarities 
with the results of the present work. Specifically, a Drude weight and coherence temperature 
scaling as 5^ are also found. The specific heat coefficient is found to scale as 1/5 in this 
case, in contrast to the present work. The reason for this difference is the existence of a 
residual entropy in the undoped spin-liquid phase of our model. However, the temperature 
dependence of the specific heat at the critical point is found to be a/T in both cases. 

Finally, we briefly discuss the possible instabilities of the metallic paramagnetic phase 
discussed in this paper. It can actually be checked that for a given J / 1, a low temperature 
and low doping regime exists in which an instability to phase separation is found, signaled 
by a negative compressibility. This is quite easily explained on physical basis for a given 
realization of the exchange couplings : the holes will tend to cluster in regions with ferromag- 
netic bonds in order to maximize kinetic energy. A proper treatment of this phase separated 
regime should take into account longer range Coulomb repulsion. In the infinite connec- 
tivity limit, an additional term VJ^ijf^i^j in the Hamiltonian reduces to a Hartree shift 
fi + V < n > of the chemical potential (thus the compressibility reads : Ky^ = i^v=o so 
that the phase separation boundary can be continuously tuned as a function of V. In future 
work, we are planning to consider other possible instabilities of this model. The issue of spin 
glass ordering [p8| does not arise for the M = oo fermionic representation considered in this 



paper |17|, but spin glass phases are indeed present for M = oo for bosonic representations 



with high enough 'spin' Even in the fermionic case, first order corrections in 1/M are 
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likely to restore a regime of spin glass ordering. Finally, an open issue is that of possible 
pairing instabilities of the metallic phase towards a superconducting state. 



B. Relevance to cuprate superconductors 

In this section, we would like to present arguments suggesting that the problem studied 
in this paper may be relevant for the understanding of some of the striking aspects of the 
normal state of cuprate superconductors. The line of arguments relies on three sets of 
experimental observations: 

• The experiments reported in in which a 61T magnetic field is used to suppress 
superconductivity strongly suggest that the ground-state of La2_2:Sra;Cu04 is actually 
an insulator, up to Sr-doping of about x ~ 0.16 corresponding to the highest Tc. This 
is true even in samples having large values of kpl, making weak localization effects 
an unlikely explanation of the logarithmic upturn of both pab and pc observed at low 
temperature. Insulating behaviour is no longer found in overdoped samples. 

• At very low doping in the La2-xSr^Cu04 compounds, a low temperature spin glass 
phase is found for x > 0.02 in agreement with theoretical arguments suggesting 
that localized holes induce locally a strong frustration in the magnetic exchange. This 
localization of the carriers induces a strong upturn of pab at low temperature in these 
samples, first in a logarithmic manner followed by an activated behaviour. Neverthe- 
less, the high-temperature behaviour of the resistivity in these samples is quite similar 
to that found close to optimal doping. 

• Inelastic neutron scattering reveals peculiar low-energy spin excitations for all un- 



derdoped samples, quite different in nature from spin waves [|l^]. For very low 
doping, these excitations occur in a remarkably low energy range, on the scale of 
lOmeV, distinctly smaller than Jaf- In a restricted range of frequency and temper- 
ature, the energy scale for these excitations is actually set by the temperature itself 



and Lu/T scaling applies |§,10| These excitations, which are present in a wide range of 
temperature (much above the freezing transition mentioned above) and in the whole 
underdoped regime, correspond to a slower spin dynamics than in a Fermi liquid, as 
is also clear from the non-Korringa behaviour of the copper NMR relaxation time. 
Similar observations have been made in the YBa2Cu306+j;Compounds |T^. This is 
particularly clear when a small amount of Zn substitution is used to suppress super- 
conductivity []T3| (we note that this simultaneously opens up again a region of glassy 



behaviour at low temperature for a rather wide range of oxygen content ||16|| ] 

In our view, these observations suggest that, in the absence of superconductivity, a T = 
metal-insulator transition occurs at some critical value of the doping x = xmi- This 
transition might be rather close to optimal doping in La2_xSr2,Cu04 0. For x > xmi, 
the incipient ground-state is a Fermi liquid, corresponding to the overdoped regime. For 
X < Xmi, the ground state is a Mott- Anderson insulator in which holes are localised at 
T = 0. At very small x (0.02 < a; < 0.05), the mechanism for this hole localization has been 



studied in Ref. |43| and involves both the freezing of hole motion due to the antiferromagnetic 
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spin background and impurity effects. This localisation induces strong frustration in the 
local exchange, in agreement with the arguments of Ref. 0. As a result, this insulator will 
have a glassy nature at T = for low doping, as indeed found in La2_a;Srj.Cu04. Beyond 
X = 0.05 however, the onset of superconductivity has prevented up to now an investigation 
of the low-temperature properties of this incipient insulating ground-state and the origin 
of the observed localization is still an open problem. It may be that the insulator looses 
its glassy character at some critical doping Xg below xmi, or that the two critical points 
actually coincide {xg = xmi)- 

As temperature is raised, the holes become gradually mobile. This quickly destroys the 
glassy ordering, leaving the system in a "slushy" state of mobile holes and spins. Neutron 
scattering and NMR experiments show that the spin dynamics in this regime is much slower 
than in a Fermi liquid state, with local spin correlations decaying (in some time range) as 
1/t (corresponding to a high density of low-energy spin excitations x"{'^) ^ const, in some 
frequency range). We view the model studied in this paper as a simplified description of 
such a slushy state of spins and holes, valid in the high-temperature quantum critical regime 
associated with the transition at T = 0, x = Xg {or xmi), as depicted schematically on 



Fig.|T2|. Indeed it is a model of a doped Mott insulator with strong frustration, in which the 
effect of quantum disordering the glassy insulating state is mimicked by taking the large-M 
limit. Fluctuations in the transverse components of the spin may actually be an essential 
ingredient in the disordering process, and this is precisely the effect which is emphasized in 
the large-M limit and produces the high density of low-energy spin excitations. 

Of course the present model is highly simplified and is meant to retain only the interplay 
of Mott localization with that of frustration in the magnetic exchange constants. As such, 
it does not include several important physical aspects of the actual materials, most notably: 

• i) The fact that frustration is a consequence of hole localisation at low temperature 
[0] (in our model frustration is introduced by hand). 

• ii) Localization of carriers by disorder (as a consequence of both i) and ii), the metal- 
insulator transition occurs at zero-doping in our model). 

• iii) The average antiferromagnetic component Ja_f of the exchange has not been in- 
cluded (in that sense we are dealing with a strong frustration limit J ^ Jaf)- This 
could be corrected for by reintroducing Jj^p in a mean-field manner, leading to : 
x{q,^^)~^ — Xioc{^^)~^ + Jaf^{(i) where A(g) is the Fourier transform of the nearest 
neighbour connectivity matrix on the lattice. We note that this formula produces a 
susceptibility peaked at the antiferromagnetic wave vector, with a correlation length 
of the order of the lattice spacing, while all the non trivial dynamics comes from local 
effects. 



For these reasons, the present model is unable to address the question of the precise nature 
of the incipient insulating ground-state of underdoped materials (or of the low-temperature 
pseudogap regime associated with it), even though the remarks made above point towards a 
phase separated regime. Various proposals have been made in the literature regarding this 
issue. One of the most widely discussed is the "stripes" picture, in which there is phase 
separation between the doped holes and the spins into domain-wall like structures. We 
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note that as long as the holes remain confined in these structures, the mechanism of Ref. 
0] implies the existence of ferromagnetic bonds in the hole-rich region, as indeed found in 
numerical calculations P3|-^. As temperature or doping is raised, a melting transition of 
the stripe structure takes place, and the model of a "spin-hole slush" introduced here may 
become relevant in the associated quantum critical regime. 

Keeping these caveats in mind, we comment on the comparison between our findings and 
some aspects of the normal state of cuprates in the regime depicted schematically on FigO: 



Low-energy coherence scale 

The present model yields a remarkable suppression of the low-energy coherence scale 
of a doped Mott insulator in the presence of frustrating exchange couplings. We find 
this scale to be of order = {5tY / J = J {5/ 5*)"^ instead of the naive (Brinkman-Rice) 
estimate 5t. We note that, with tj J 5, and J ~ 12007^, this scale can be as low as a 
few hundred degrees. If relevant for cuprates, this observation suggests that the normal 
state properties may well be associated, over an extended (high-) temperature regime, 
with incoherent behaviour characteristic of a quantum critical regime dominated by 
thermal effects. We note however that the present model, as any model in which low- 
energy excitations are local in character, would lead to a large effective mass at low 
temperature, directly proportional to In cuprates, additional physics sets in at 

lower temperature (cf iii) above) which quench the corresponding entropy, leading to 
the experimentally observed moderate effective mass . 



Photoemission 

In the incoherent "slushy" regime T > e*p^ we find a single particle Green's function 
decaying as 1/\/oj (and an associated single-particle lifetime ImS oc -Joj), leading to 
a markedly non-Fermi liquid tail of the photoemission intensity. It is worth noting 
that precisely this form has been recently shown to provide a rather good fit to the 
high-frequency part of the photoemission lineshape above the pseudogap temperature 
in underdoped Bi2Sr2CaCu20s+x HH- It has been recently argued that the I/a/cJ 
behaviour also holds in the strong coupling limit of antiferromagnetic spin fiuctuation 
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Resistivity and Optical conductivity 

We would like to emphasize again the mechanism which yields a linear resistivity in the 
incoherent regime of our model, starting from a single-particle self energy behaving as 
^/uJ. This holds when scattering is local and incoherent so that the effective quasipar- 
ticle bandwith (dispersion) can be neglected in comparison to lifetime effects. In this 
limit, conductivity should be thought of in real space as a tunneling process between 
neighbouring lattice sites. This mechanism has a higher degree of generality than the 
specific model considered in this paper, and should also apply to other models in which 
the same power-law behaviour of the self-energy holds, such as the model of Ref. 
p5|] . This model has been proposed in connection with the normal state properties of 



underdoped cuprates above the pseudo-gap temperature . 

We note that the magnitude of the linear resistivity in this incoherent regime is larger 
or comparable to the Mott limit (a/i/e^), as is actually the case over a rather extended 
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high temperature regime in underdoped La2-x^^xCu04, and is a quite general 
feature of 'bad metals'. 

Regarding optical conductivity, the form we have obtained is quite similar to the 
Marginal Fermi Liquid one, which has been shown to provide a very good fit to 
the data of e.g. Ref. 



Neutron scattering 

Neutron scattering experiments on non-superconducting Lai.95Bao.o5Cu04 and 
La2_a;Srj.Cu04 with x = 0.04 ||T0|JT3[| have revealed spin excitations which are centered 



at the wavevector Q = (vr, vr) with a rather large momentum width. The frequency 
dependence of these excitations display co/T scaling and have been successfully fitted 
by scaling forms very similar to the one found in the present model |]8|JT0|. At higher 
Sr concentration, one of the most notable feature of the neutron scattering results is 
the appearance of sharp peaks at incommensurate wavevectors. It is likely however 
that these peaks only carry a small fraction of the total spin fluctuation intensity, as 
suggested in particular by comparison to NMR data. A broad, weakly g-dependent 
contribution most probably persists up to high temperature, carrying a large part of 
the total weight, and hard to distinguish from "background" noise in neutron exper- 
iments [|^. In YBa2Cu306+y, suppression of superconductivity by Zn doping allow 
to investigate the spin dynamics of the normal state down to low temperature |l5| , p!4 



Apart from a very low temperature quasi-elastic peak (associated with spin freezing 
into spin-glass like order), neutrons scattering results for y = 0.39 reveal a strong en- 
hancement of low-frequency spin fluctuations at low temperature, with a distinctively 
low energy scale and a strong temperature dependence down to very low temperature 
(compatible with u/T scaling in a limited range). These features are qualitatively 
similar to the low-energy excitations found in the present model. There is furthermore 
experimental evidence [|15[ that these low-energy excitations are associated with the 
disordering of the spins by transverse fluctuations, as in our model. 
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APPENDIX A: DERIVATION OF THE SADDLE-POINT EQUATIONS 

In this appendix, some further details on the derivation of the saddle-point equations 
in the large-M limit for the single site model defined by Eqs.(|lT],|12D are provided. In the 
following, we will drop the index a in G"^*(r - r') = -jj < ifab^'')ir){pJ'b^){T') >. In (|g) 



24 



the brackets denote the average with the action specified in subscript. As the action S is 
invariant under translations in imaginary time and under the action of SU{M) (the rotation 
invariance for M = 2), C and R take the following form : 

K^sir, r') = -d^sSp^R^'ir - r') + d^pd^sR^'ir - r') (Al) 

The quartic term in / in (^) is decoupled using a bi-local field P"-^{t,t'). Using the 
expression of the spin operator Eq. (H) and the change of variable : 

6(r) = ^0(r) (A2) 
the single-site partition function can be rewritten as : 

^sin.,e-site = / V<P^ Vcj^V XV P 6- (A3) 

with the actions : 

•^1 = 2/^^^ ^^"i^)drrir) - log ^0 + Y E // drdr' R-\r - t')P-\t, T')P'\r\ r) 

1 ab 

- J2 



S\ = \Y.[jd'^{^-'P\'^)m)) (A4) 
In this expression, Zq is defined by 

Zo[(/),P,A] = j p/tp/e-^«°['^-P'^'^l (A5a) 



with : 



a •' 



^'"tW(t) - 1 



2 / 



a,b 



E / / drdr' /t'^(r)0"(r)G"(r - r')0^"(r')r(r') (A5b) 



In the limit M — ^ oo, Zgingi<,_sito is controlled by a saddle point with respect to P°'^(r', r), A(r) 
and 0(t). We assume a condensation of the boson : after the change of variable ( |A^ ) is 
taken to be a finite constant at the saddle-point : 0sp(t) = and A is static : iXgpir) = Xq. 
Moreover, in this limit the correlation functions of / are given by the average with the action 
Sqo taken for these values of P, A, 0. As 5*00 is quadratic in / (and the boson is condensed), 
the model is completely solved in this limit as soon as Gf has been calculated. The saddle 
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point equations are given by the minimisation of 5*1 with respect to P''°(r', r), A(r) and 0(r) 
respectively, which leads to : 

l = 6-2<nT)p%r) >Soo 



Jo 

and finally gives Eg. ( |14a| , p.4b| , |14c|Jl4d| ) given in the text. 



(A6) 



APPENDIX B: LUTTINGER THEOREM 



In order to find the volume of the Fermi surface in the interacting system, we proceed 



along the lines of Ref. ||32| and we observe that the auxiliary fermion self-energy can be 
obtained as the functional derivative of the following functional: 

5$ 



^ = j dt {G'f{t)G'f{-t)f , Sf (t) 
The number of particles reads: 



6Gf{ 



-t) 



-oo 2m 



G'fiu) e 



and we use the identity: 
d_ 



/oo 
dtDit) ln(cj + /i 
-oo ^ 



Ar 



(Bl) 



(B2) 



(B3) 



Using the invariance of the Luttinger-Ward functional under a shift of all frequencies 
{GiiS) Giio + i^)), the integral of the last term vanishes: 



dio 







(B4) 



and the integral of the first term can be explicitly calculated by transforming to retarded 
Green's functions (denoted by ) in the following manner: 



1-5 



deD[e) 



duj 



— 6'^lnG^(e,a;)e 



-00 2ni 



'Gfie, 



Gf{e,u;) 



(B5) 



As has no pole nor zeros in the upper half plane, the first integral can be closed there 
and vanishes. Hence we have : 
1-5 i-c-o 



deD{e)Q (fi-Xo- S/(zO+) - 5e 



From Eq. ([B6|) and the definition of yUo, we finally obtain: 

/i(T = 0) - Ao(T = 0) - = 0, T = 0) = 5fXo{5) 



(B6) 



(B7) 



which is the desired relation and insures that Luttinger theorem holds in the presence of 
both the constraint and the magnetic scattering. We also checked that this property is 
verified in our numerical calculations at T = 0. 
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APPENDIX C: NUMERICAL METHOD 



In this appendix, we explain the main steps that we followed in solving numerically the 
saddle point equations Eqs.(|T3). 



1. Computation of the Green function Gf{uj,T) 

The calculation of the Green's function is divided in two steps. First a Matsubara 
frequency/imaginary time algorithm is used in an iterative manner in order to find the value 
of the chemical potential fi and Lagrange multiplier A for a given doping 6, interaction 
strength J/t and temperature. Convolutions are calculated using a FFT algorithm and a 
simple iteration is used: starting from a given G, a self-energy is obtained which is then 
reinjected into the expression for G until a converged set (G, S) is reached (for given values 
of fj,, X). A second routine uses the previous one to adjust /x and Aq in order ( p.4c| , |14d| ) to be 
satisfied. 

Once the imaginary-time Green's function and values of fi, A have been found using this 
imaginary time algorithm, a different algorithm is used to obtain real frequency Green's 
functions and spectral densities. This is done in the following manner. We consider a 
finite-temperature generalisation of the Green's function with the Feynman prescription: 

G{uj) = {l-nF{u))G^{uj)+nF{uj)G^ (CI) 

which reduces to the usual Green's function Gp ^.t T = 0. We now define: 

S(t) = J^G{tfG{-t) (C2) 

where t is the real time. Expressing both S and S'^ as integrals of the spectral density with 
the spectral representation of G^, we obtain, after some calculations (R is a superscript 
denoting retarded quantities) : 

= ^^(uj) + 2z7rJ^ JJ duidu2np{ui)p{uJi)np{ijj2)p{ijJ2)nF{'^ — cui — uj2)p{uj — uji — 0J2) 

(C3) 

where p = ImG^ and np is the Fermi factor. 

At T = 0, Eq. ( P3| ) shows that S simply coincides with E^, the usual self-energy at 
T = 0, with the Feynman prescription and that Eq. (|I^) can be rewritten as, at T = : 



{G}\uj))-' = uj + pi-Xo~ (tSfG'^iuj) - ' 



mt) = J\G}it)yG}i-t) (C4) 

together with the equation corresponding to (|14c| , |14d| ). Note the change of the sign in 
front of J^. This form was used in (pSf). From ( |C4| ) one can write an algorithm for the 
computation of G^ in the T = formalism, similar to the computation in imaginary time. 
Of course, as our large-M limit performs a resummation of the perturbation theory, one can 
also obtain these equations with the diagrammatic rules, but these rules do not apply at 
finite temperature to the Green's function G in a systematic manner. 
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A finite temperature, we use the following iterative algorithm : Starting from G^, we 
get G. We then obtain E by direct convolution in real time and the second term of the r.h.s 
of ( |C^ ) by a double convolution of npp. Hence we obtain S-^ and go back to with Eq. 

(0). 

As a starting point of the iteration, in order to speed up convergence, we take an analytic 
continuation of the solution in imaginary time, obtained by a standard Pade approximation. 
Note that the parameter and Aq are fixed in this iteration on the real axis, since they have 
been calculated before in the Matsubara formalism. 

As soon as the Green function has been obtained, some other quantities are straightfor- 
wardly calculated from the formula given in the text. In particular, x"^ is expressed as a 
convolution. 

However the computation of the uniform susceptibility x (considered in Appendix ^ is 
more involved : we solve Eq . (|E3|) for g by another iterative loop analogous to the previous 
ones. 

The scaling function describing the effect of the doping at T = 0, is computed from (pH]) 
by an iterative algorithm similar to the previous ones. 



2. Computation of the resistivity 

We also give some useful details on the numerical calculation of the frequency-dependent 
resistivity. There, it is very convenient to integrate analytically over e in (^) using (^0[) and 
the relation : 



de 



Die) 



GfM 



Aiu)-e ■ 

We thus obtain (the Green function is the retarded one) : 



(05) 



Re aiuj) 



6^ 



duKe 



G Av + uj) - G fiv) G Av) - G fiv + uo) 



For the dc-conductivity, Eq . (|C^ ) simplifies to : 



_ A{uj + z/) - A{u) A{v) - A{u + uj) 



(C6) 



a(T) 



327r2T 



dio 



Re 



b'^G){u) 



cosh^ f 



(C7) 



In both case the integrals are computed simply by transforming them into a Riemann sum. 



APPENDIX D: SCALING ANALYSIS 

In this Appendix, some details about the thermal scaling analysis of Section |III C| in the 
spinfluid regime are provided. As explained in the text, in this regime e*p disappears from 
the thermal scaling functions and thus the calculation can be performed in the undoped 
model (5 = 0. In this Appendix, Gf and will denote the thermal scaling function of these 
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quantities. They satisfy the scaled saddle point equation G f{u!/T)~^ = —Tif^u/T). The 
calculation is very similar to the low temperature, low frequency analysis made in so we 
here just give the main steps of the analysis. We first check that the scaling behaviour ( |50[ ) 
in imaginary time solves the saddle point equation (for 6 = 0) using the Fourier formulas 



(which follows from 3.631) : 



Gf{iuj„ 



3(jT)-i(-i)"r( 



^n/3 \ p 3 _|_ UJnl3 

27r / V4 27r 



(Dla) 



miv^(-lfr(-i) 



^nl^ \ p / 1 _|_ ^n/3 \ 



(Dlb) 



where Un = {2n + 1)ttT are the Matsubara frequencies. One can then show that ( ^91) and 
(pT|) are the scaling function on the real axis using the following method. First we use 



GAr) 



+00 



—^Pf{e)de 0<r</? 



:d2) 



and 



+00 



dt 



TX 



-itu 



cosh(7rt) 



(27r 



A £U 

2 27r 



r(A) 



< A < 1 
u real 



(D3) 



(See formula 3.313.2 of We then deduce the full Green function (and thus S/) by 

performing the Hilbert transform of p/ using ( P3| ) again and 



dx- 



(sinh 

(See formula 3.112.1 of [H 



-i/? r(f-g)r(i -A) jo<A<i 
- r(i-f-z( 



z real 



(D4) 



APPENDIX E: UNIFORM SUSCEPTIBILITY 

In this appendix, we briefly explain how to calculate the uniform susceptibility x and 
analyse its behaviour at small temperature in the undoped model 5 = 0. 



1. Effect of a magnetic field 

The magnetic field is introduced in the SU (M) Hamiltonian (^ in the following way : 

6H = -h [flh - /2V2) (El) 

This formula clearly reduces to the usual one for M = 2. Here, only color 1 and 2 are coupled 
to this field but note that this choice is not unique though convenient for our calculation 
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(more generally the magnetic field must be coupled to an element of a Cartan subalgebra of 
SU(M)). The large d and large M limit computation is similar to the zero field one explained 
previously, although it is more involved. A simplification occurs in this double limit : due to 
the fact that only 2 colors over M are coupled to h, the Green function for colors i > 2 
are solutions of the zero-field equations (p^. Moreover we find = G and = G^^ 
where G^ is given by : 

G'}{iuJn)) ~^ = ioOn + 1^- Xo + h- {t5fG){iUn) - T.){iujn) (E2a) 



S5(r) = -J^G){t)G)=\t)G)=\-t) (E2b) 

(/i, Ao and G^"^ are always determined by Eqs. ([T^ ). The magnetisation is given by m = 
G}(0-) -Gf(O-). Let us define g by G)-G) = hg + 0{h^). From Eq. ( |E2i| , [E2b|) , g satisfies 

{tdy - Gj [tLUn) 

Kir) = 25{t) + .Pg{T)Gf{r)Gf{-T) (E4) 
With these notations, the uniform susceptibility is given by % = g{T = 0^). 



2. Low temperature behaviour of x for the undoped model 

In the undoped case, Eq.(^3|) reduces to g{uj) = —G'j{uj)K{uj) and the susceptibility at 
T = is formally given by : 

X(T = 0) = r g{u;)du; (E5) 

From the low frequency behaviour Eq . (|T7|) we have g{u}) ~ —i const. K{uj)/uj. Thus we 
have to investigate the low-frequency behaviour of K. Generalising the Luttinger theorem 
(as expressed by Eq.(0)) for the colors 1 and 2, we obtain (at T = 0) : 

fi±h- Xo-J:y\uj = 0) =6fi^\6) (E6) 

In this equation, fJ^o'^iS) is given by : 

deD{e)=n^ (E7) 
where is the number of particles of color 1. Hence we obtain : 

2 - (4(0) - Sj(0)) = 5/M_/fi!(^ (E8) 

Taking first the limit 5^0 and then h —>■ Owe obtain finally (as fiQ is bounded by definition): 

K{u = 0) = (E9) 

Thus the leading low- frequency singularity in g cancels so from Eq.(^) we see that x is 
smaller than InT at small temperature. Strictly speaking we can not prove from the previous 
argument that x reaches a finite value at zero temperature, but it is a very natural guess 



which is moreover very well supported by our numerical calculation as displayed in Fig.|13 
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FIGURES 
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FIG. 1. Crossover diagram as a function of temperature and doping. The coherence scale (*p 
is indicated by a dashed hne and is given by e*p ~ J{5/5*Y for 5 < 5* , e*p ^ 5t ior 5 > 5* , with 
5* = J/t. Below e*p, Fermi liquid behaviour holds. For 5 < 6* , an intermediate "quantum critical" 
regime is found in the range e*p < T < J, in which charge transport is incoherent and spins have 
a marginal Fermi liquid dynamics. 
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FIG. 2. Local dynamical susceptibility x(^^(a;, T = 0) of the undoped spin liquid. Inset : 
spectral function. 
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FIG. 4. T = scaling function associated with the spectral density Pf{io) = j§4'fi'^/^*F) 
the low-doping regime. 
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FIG. 5. Physical electron quasi-particle residue vs. doping for J = 0.3, 1,2.5 (the propor- 
tionality factor 2/M has been set equal to 1) 




FIG. 6. Conduction electron spectral density p{e]^, io) for 6 = 0.04, T/t = 1/300 and J/t = 0.3 
and for 2 values of the energy ej^. The arrow indicates the cross-over between the Fermi-liquid 
regime and the spin-liquid regime, as explained in the text. 
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FIG. 7. Spectral functions Trpf{uj) for 6 = 0.04 and J/t = 0.3 (correspond- 

ing to Gp/J — 1.8 10"^). The different curves correspond from top to bottom to 
T/t = 1/200, 1/50, 1/25, 1/10, 1. Inset : thermal scaling function Eq. 
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FIG. 8. Local dynamical susceptibility Xiod^) S = 0.04 and J/t = 0.3. The different 
curves correspond from top to bottom to T/t = 1/200,1/50,1/25,1/10,1/5,1. In the temper- 
ature range e*p < T < J and for frequencies to < J, these curves scale on the universal form 
Xro>,T) = ^tanh# 
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FIG. 11. Real part of the optical conductivity Rea{uj) versus u, for 6 = 0.04 and J/t = 0.3. 
The different curves correspond to T/t = 1/200,1/100,1/50,1/25,1/10. Here e*p/t = 1.8 lO^^. 
Inset : the curve corresponding to T/t = 1/100, plotted in log-log coordinates, in the frequency 
range T ~ < u; < J. The 1/lo behaviour described in the text is clearly visible. 
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FIG. 12. Schematic crossover diagram for cuprates illustrating : i) the existence of a 
metal-insulator transition as a function of doping at T = and ii) the possible relevance of our 
model to the corresponding quantum critical regime (other features like the pseudogap and the 
Neel temperature have not been depicted) 
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FIG. 13. l/xiT) (solid line) and l/xioc{T) l/InT (dashed line) versus T as calculated 
numerically in the undoped model. 
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